Size Dependence In The Disordered Kondo Problem 



Ivar Martin, Yi Wan, and Philip Phillips 

Loomis Laboratory of Physics 
University of Illinois at Urbana- Champaign 
1100 W.Green St., Urbana, IL, 61801-3080 

Abstract 

We study here the role randomly-placed non-magnetic scatterers play on 

the Kondo effect. We show that spin relaxation effects (with time t°) in the 

vertex corrections to the Kondo self-energy lead to an exact cancellation of 

the singular temperature dependence arising from the diffusion poles. For a 

thin film of thickness L and a mean-free path £, disorder provides a correction 

to the Kondo resistivity of the form t° /(kpLi 2 ) InT that explains both the 

disorder and sample-size depression of the Kondo effect observed by Blachly 

and Giordano (PRB 51, 12537 (1995)). 
PACS numbers:72.10.Fk, 72.15.Nj, 75.20.Hr 
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At low temperatures, the resistivity of a metal alloy acquires a logarithmic temperature 
dependence [|1| in response to spin-flip scattering between local magnetic impurities and the 
conduction electrons. This behaviour persists down to a temperature (the Kondo temper- 
ature, Tfe) where the magnetic impurities and conduction electrons begin to condense into 
singlet states. While perturbation theory is sufficient to establish the existence of the InT 
term, its presence ultimately signifies that perturbation theory is breaking down. Alterna- 
tively, spin-flip scattering between conduction electrons and localized magnetic centers has 
a singular frequency (u) dependence. Magnetic impurities are not alone in this respect. 
It is well-known that even non-magnetic impurities can generate a singular (In a; in d=2) 
frequency dependence in the conductivity 0. In a sample containing both magnetic and 
non-magnetic impurities, the question arises: which singularity ultimately wins or can the 
interplay between the singularities lead to a suppression of either localization or the Kondo 
effect? In this letter, we resolve these questions. 

The motivation for this study is two- fold. First, while there have been numerous treat- 
ments of this problem ||- |7|], a clear consensus has not been reached. To illustrate, Everts 
and Keller || were first to show that non- magnetic scattering contributes a 1/VT to the 
Kondo self-energy that dominates the Kondo InT singularity at low temperatures in d=3. 
Bohnen and Fisher || argued, however, that such a term would not survive in the con- 
ductivity. More recently, Ohkawa, Fukuyama, and Yosida || showed that disorder results 
in a singularity of the form T d l 2 ~ 2 in the conductivity. At low temperatures this singular- 
ity dominates the Kondo InT. As a result, these groups conclude that static disorder can 
mask the Kondo resistivity as T — > 0. On the experimental side, Blachly and Giordano || 
recently measured the conductivity in a series of thin films containing magnetic as well as 
non-magnetic impurities. They found no evidence for the T d l 2 ~ 2 singularity but observed 
instead a suppression of the Kondo resistivity as the strength of the disorder increased. 
Earlier experiments by Korn || also failed to observe the T d ^ 2 ~ 2 singularity but observed 
instead an enhancement in the Kondo resistivity. The point of agreement between these 
experiments is that disorder couples non-trivially to the Kondo effect and ultimately modi- 
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fies the coefficient of the InT dependence. Given the strong dimensional dependence of the 
localization transition, disorder could eventually lead to a sample size dependence of the 
Kondo effect. 

At the outset, we set aside the still controversial issue (ref. 8c) of the sample size 
dependence and focus on the seemingly straightforward problem of the role non-magnetic 
disorder plays in the Kondo effect. The new wrinkle we introduce in this problem is the 
feedback effect spin scattering has on localization. While it is standard to consider the 
direct influence of localization on the Kondo effect, the reverse effect has not been included 
|T0| . Nonetheless, it is well-known that electron scattering by disordered Heisenberg spins 
introduces a cutoff of the diffusion pole in both the particle-hole (diffuson) and particle- 
particle (Cooperon) channels except for the S = particle-hole channel When spin- 



scattering is included in the diffusion propagators, the fate of the T d ^ 2 ~ 2 singularity rests on 
whether the S = particle-hole propagator contributes to the Kondo self-energy. We show 
explicitly it does not. 

The starting point for our analysis is a model Hamiltonian H = H + H s d that contains 
both normal impurities 

Ho = J2 (e* - e F ) ala ka + 1 £ e^'^ala^ (1) 

ka k,k',i 

as well as magnetic scatterers 

H sd = ~ J2 e l{k - k ' >Rn ■ S n aU*v. ( 2 ) 

R n ,k,k',(T,a' 

The operator a\ a creates an electron in a plane wave state with momentum k and spin a and 
energy v measures the strength of the scattering with the non-magnetic disorder, J is 
the exchange interaction, R n denotes the position of the impurities, magnetic or otherwise, 
S n is the spin operator for the magnetic impurity at site n, a is the Pauli spin operator, and 
Q is the volume. The two natural timescales in this problem are, r° and r , the magnetic 
and non-magnetic scattering times. In terms of the density of states of the host metal, p Q 
and the concentrations of magnetic and non-magnetic scatterers, n s and n Q , respectively, 



we have that h/2r° = 37m s p \J\ 2 /4 and h/2r = 7m p \v\ 2 . The total scattering rate 
is 1/t = l/r° + l/r G . To measure the strength of the non-magnetic disorder, we define 
A = h/(27TSFT ). We assume that the concentration of localized spins is dilute so that long- 
range spin glass effects are irrelevant. Also, we work in the regime in which normal impurity 
scattering dominates, l/r G ^> l/r° 

To evaluate the conductivity above we must first calculate the Kondo self energy. To 
include the dynamical effects of the localized spins, it is sufficient to calculate the self energy 
to third order in the exchange interaction J. At this order, static disorder can be included by 
decorating the single and double spin-flip vertices with Cooperon and diffuson propagators 
0, 0. In previous work ||, |7|], spin- independent Cooperons and diffusons of the form 
C(Q,u>) = D(Q,u) oc ( D q2_.^) were used where Q and u are the net momentum and energy 
transfer and D = 2heFr/dm is the diffusion constant. However, this is inconsistent because 
C(Q,u) and D(Q,u>) are coupled to electron lines of different spin. Such propagators are 
well-known [[□]] to depend on spin and hence we include explicitly the spin dependence 
here. If all scattering processes are treated in the first Born approximation, the Cooperon 



propagator JTTJ] is transformed to 



h 2 

C a f3~f8 = ~ o7n7)2 : , o/ o\ (<W<V _ a <*P ' °V) 

8np T 2 (DQ z — iuj + 2/t°) 

h 2 

+ 8n Po r 2 (DQ 2 - i u J + 2/3T°) + ' (3) 



and the diffuson becomes 



h 2 



Da/3-yS - ~ 2( n ^, 2 (£a/3(V + Va(3 ' Vys) 

8np T z (DQ 2 — VjJ) 



h 2 



87rp T 2 (DQ 2 - iuj + 4/3r°) 7 p 7 



(4) 



The survival of the diffusion pole in the spin- dependent diffuson is a consequence of particle- 
hole conservation. The two terms in each of these propagators correspond to singlet and 
triplet scattering, respectively. Pairs of spin indices a/3 and are indexed chronologically 
along the particle lines that comprise the diffusion ladder diagrams. 



The diagrams shown in Fig. la contain the dominant quantum corrections to the Kondo 
self-energy at third order in the presence of disorder. The sum of all such diagrams is 

2 

^3 q (k, ie n ) = — J2 V a/3vr)( iuJ ei iuj m )G(ie n + iu m , q) x 

G{ie n + iuj£, k + Q) + n \v\ 2 '^2 G 2 (ie n , k')G{ie n + iu^, k' + Q) x 

V 

{D aa ^{iu) £ , Q)D luria (iu e , Q) + C aajv (iu £ , Q)C/3 irja (iui, Q)) (5) 



where G(ie, q) is the electron Green function 



G(ie,q) 
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(6) 



ie + ep — h 2 q 2 /2m + i(h/2r)sgn(e) ' 

the electron energies are the Matsubara frequencies, e n = (2n + 1)tcT, the psuedofermion 
energies are = (2k + 1)tcT, uj^ = 21ttT, DQ 2 < h/r and (e n + lu/)U£ < 0. We have set 
ks — 1- The factor of 2 arises from the two possible couplings of the diffusion progagators to 
the internal electron lines and the ± from the two orientations of the psuedofermion loops. 
The psuedofermion part involves a trace over the components of the impurity spin operators 
and hence simplifies to 



;J 3 n s p 



- — {Smo — <W)(1 ~ $io) + - — Sio(l — 5 m o) ± — 5 m o5io 
uoi iu) m 2 



I a a \ 
\°af3 u vr)) 
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From the psuedofermion contribution, we see that the sum over the spin indices separates 
into two identical sums of the form, D^a^a^. If we use the identity Y^ap^va^p-y)'^ 
— cr" 7 , we find immediately that the cancellation of the divergent diffusion terms 



a 
'a/3 



a/3 a/3 
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from the 3 rd order Kondo self-energy is exact. To any order in J, the cancellation of the 
diffusion pole can be seen as follows. In the most divergent approximation, each diffuson 



encircles a vertex that is exactly equal to the Abrikosov |lz] vertex function r oc a- S. When 
this function is now multiplied by D s=0 and summed over the spin indices, the cancellation 
to all orders follows immediately from Eq. (H). This is one of the principal results of this 



paper. The cancellation of the S = component of the diffuson is fundamentally tied to the 
fact that the Kondo interaction does not conserve spin. Further, it signifies that the resultant 
conductivity is independent of the pure charge density propagator. Summing over the spin 
indices in the remaining propagators in the self energy reduces the problem to one in which 
the diffuson and Cooperon are spin independent: D = h 2 / (2np T 2 )(DQ 2 - iuo + 4/Srf)- 1 
and C = h 2 /(47tp T 2 )[(DQ 2 - iuj + 2/r°)- 1 + (DQ 2 - iu + 2/3r s °)- 1 ]. When l/r s ° = (or 
equivalently, T ^> % / t°) , we recover the standard form for these propagators. 

To calculate the resistivity, we evaluate the standard self-energy as well as the Cooperon 
weak-localization diagrams || . Because the results of the calculation are rather lengthy, we 
present here only the asymptotic behaviour. In the limit T ^> fi/r° (as in the case when the 
magnetic impurities are dilute), we recover the inverse temperature dependence 

h h —7ihp \j h , T ft 

< - Po XJ- (9) 



of refs. ( 0, |)|). In the opposite regime, T ^/t", the limiting forms of the Cooperon 
and diffuson relaxation times 

h h (h 31n3\ . r h . h 

+ 7^a = ~ ( o + — I PoAJ- In — (10) 



2t d 2t c \2 4 J ro t Tr b 

are both logarithmic functions of temperature. 

The final contribution to the relaxation time comes from the the Cooperon weak- 
localization diagram. In two dimensions in the presence of spin-flip scattering, the weak- 
localization contribution is Acr/ oc = —e 2 /(27c 2 h) ln(\/3T es /T ), where h/2r es = 8h/(3r°)(l — 
p Jln(ep/T)). Physically, t €S plays the role of the inelastic scattering time in the weak- 
localization correction. Inclusion of the 3 rd order correction to the spin scattering time, 
enhances the spin-flip scattering rate, thereby weakening the effects of localization. To see 
this more clearly, we expand the argument of the logarithm for temperatures well above the 
Kondo temperature: 



We see clearly that the Kondo interaction reduces the weak localization correction because 
J < 0. 

We collect all the contributions discussed above to determine the conductivity. In the 
temperature range C T < h/r°, Cooperon, diffuson, and weak-localization corrections 
are logarithmic in temperature. Combining the results from Eq. (|TUD with the weak- 



localization correction, we find that the magnitude of the logarithmic part of the conductivity 



r 4r p J / a\ t s\ ,_ £f 



Aa+=a ^(l + lAXfJlnf (12) 

is enhanced by disorder. The first term in this expression arises from the unperturbed 
Kondo effect and the latter from the interplay with disorder. Inclusion of disorder in the self 
energy always enhances the Kondo resistivity by increasing repetitive scattering at magnetic 
impurities. 

For temperatures T ^> h/r°, the self-energy contribution to the relaxation time scales 
as 1/T, whereas the weak- localization correction is proportional to InT. However, com- 
parison of the magnitude of these corrections (see Eqs. (^) and (|TT])) reveals that the 
weak-localization term dominates and the magnitude of the resultant logarithmic correction 

A^ = a ^(l-^)m- (13) 

is suppressed by the disorder. The ratio A/r Q scales as l/£ 2 , where i is the mean-free path. 
We see then that in the dilute impurity regime, disorder suppresses the Kondo effect. The 
crossover from enhancement to suppression of the Kondo effect occurs because the magnitude 
and functional dependence of the quantum corrections to the self-energy are determined by 



the shortest of two length scales: the phase-breaking length, = JDr°/h and the diffusion 



length, Lt = JD/T. The latter arises because coupling of diffusion propagators to internal 
electron lines in the self-energy leads to an effective electron-electron interaction. 

Let us now apply our results to the experiments on thin films by Blachly and Giordano 
||. In all of their samples the film thickness L satisfied the inequality i < L L^. Hence, 
we can treat the films as quasi-2d with respect to localization, but because I < L the electron 



gas is characterized by a 3- dimensional density of states p Q = l/(27r) 2 (2m//t 2 ) 3 / 2 e 1 / 2 with 
a diffusion constant given by D = 2TiEpT /?>m. The summation on Q in the Cooperon 
and diffuson is restricted to small momentum transfers such that DQ 2 < l/r . However, 
for thicknesses of the sample on the order of £, the smallest wave vector in the transverse 



direction does not satisfy this constraint. To rectify this problem, Volkov |T3[ showed that 
surface boundary conditions must be treated consistently. For thin films, his treatment 
shows that the boundaries always give rise to a strictly 2-dimensional weak-localization 
correction and an explicit finite size dependence. To account for the former, the momentum 
integration in the Cooperon and diffuson must be restricted to the plane. The density of 
states that arises from converting the sum to an integral will be the 2-dimensional density 
of states p 2D = Tcp / (kpL). Hence, the self-energy diagrams will generate a size-dependence 
to the conductivity. The explicit finite-size weak-localization correction is |13[ Aer/ oc = 
-e 2 /(2n 2 hL) In (VSt €S /t (smh(L/£)(£/L))). The size-dependence in the logarithm yields 
an effective size dependence in the spin-relaxation time. However, this will not affect the 
temperature dependence of the conductivity. Hence, the only size dependence that is coupled 
to the temperature is the 1/L prefactor of the weak- localization correction. 

We now combine these results in the low and high-temperature limits discussed earlier. 
In the two limits, we obtain 

Aa r = | (l + ln ¥ if T k «T< h/T° 

1 °o*^ (l - ^tp) Inf if T k ,h/r° « T 

an explicit size and disorder correction that scales as 1/(£ 2 L). In the concentrated impurity 
limit T < h/r°, increasing disorder enhances the resistivity. In Cu(Fe) alloys at impu- 
rity concentrations ranging from 0.3 — 2.1%, Korn observed an enhancement in the Kondo 
resistivity that is consistent with the first equation above. However, in the dilute limit, 
T ^> h/r°, we predict a suppression of the Kondo effect as the disorder is increased and the 
size of the sample decreases. In the experiments of Blachly and Giordano ||, h/r° m 0.1K 
which is much less than the Kondo temperature for Cu(Fe). The second of equations should 
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be valid. Fig. 2 shows a comparison between the experimental data and the theoretical pre- 
dictions. The best fit to the data was obtained with r° = 0.52ns which is consistent with the 
experimental range of 10~ 10 s. As is evident, theory and experiment are in good agreement. 
We also obtained quantitative agreement with the experimental data when the sample size 
was varied. We conclude that disorder can suppress the Kondo resistivity and give rise to a 
sample size dependence of the form 1/(£ 2 L). We note in closing that a recent theory of the 
Kondo size dependence in clean samples has been proposed by Ujsaghy, Zawadowski, and 
Gyorffy [Q. This approach applies strictly in the ballistic case where £ 3> L. 
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FIGURES 

FIG. 1. Feynman diagrams contributing to the Kondo self-energy (£). The dashed lines 
correspond to Abrikosov psuedofermions and the double solid lines to diffusons and double dashed 
lines to the Cooperons. The Greek letters indicate the spin. The X indicates a single non-magnetic 
impurity scattering event. 

FIG. 2. Comparison of the theroretical prediction for the Kondo resistivity predicted from the 
second of Eq. (|lj|) with the experimental data of Blachly and Giordano ( ||) Fig. 7. The horizontal 
axis measures the strength of the static disorder through the mean-free path. 
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Fig. 1a. 



